Quasiattractors in coupled maps and coupled dielectric cavities 
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We study the origin of attracting phenomena in the ray dynamics of coupled optical microcavities. 
To this end we investigate a combined map that is composed of standard and linear map, and a 
selection rule that defines when which map has to be used. We find that this system shows attracting 
dynamics, leading exactly to a quasiattractor, due to collapse of phase space. For coupled dielectric 
disks, we derive the corresponding mapping based on a ray model with deterministic selection rule 
and study the quasiattractor obtained from it. We also discuss a generalized Poincare surface of 
section at dielectric interfaces. 

PACS numbers: 05.45.-a, 42.55.Sa 



I. INTRODUCTION 

An attractor in classical dynamics is, colloquially 
speaking, a (invariant) set of points to which all trajecto- 
ries starting in its neighborhood (more precisely, its basin 
of attraction) converge. More precisely, an attractor can 
be defined as a closed set A with the following properties 
PHI. First, A is an invariant set: any trajectory x(t) 
that starts in A stays in A for all time. Secondly, A at- 
tracts an open set of initial conditions: there is an open 
set U containing A such that if x(0) G U, then the dis- 
tance from x(t) to A tends to zero as t — )• oo. This means 
that A attracts all trajectories that start sufficiently close 
to it. The largest such U is called the basin of attraction 

of A m. 

The fact that volumes have to be conserved in con- 
servative systems implies immediately that they display 
no attracting regions in phase space [5|, Q. However, a 
quasidissipative property has been reported in a com- 
bined map, namely a piecewise smooth area-preserving 
map which models an electronic relaxation oscillator with 
over- voltage protection [7[. The quasidissipative prop- 
erty eventually converts initial sets into attracting sets, 
which was correspondingly referred to as the formation of 
a quasiattractor. As a result, this system behaves partly 
dissipative (outside the quasiattractor, i.e. before trajec- 
tories have reached the quasiattractor) and partly con- 
servative (inside the quasiattractor). To the best of our 
knowledge, the existence of quasiattractors has been ob- 
served so far only in the piecewise smooth area-preserving 
map. In this paper we show that coupled dielectric cavi- 
ties are another system class with this property. 

Each optical interface is characterized by the splitting 
of an incoming ray into transmitted and reflected rays 
which is repeated at each reflection point. For optical 
systems consisting of more than one dielectric building 



block such as the combination of two disks, the ray split- 
ting dynamics is especially complicated. A ray model 
with deterministic selection rule (RMDS) was proposed 
J8j to effectively describe the resulting dynamics. Its jus- 
tification is underlined by a nice agreement with wave 
calculations [8[. 



A noticeable and characteristic property of ray dynam- 
ics in the RMDS is that all initial rays eventually arrive 
in a certain region of phase space - namely an island 
structure |8J|. In other words, attractors occur. Their ex- 
istence is, on the one hand, a surprise because the under- 
lying ray dynamics is Hamiltonian. On the other hand, 
wave solutions for the individual resonant modes are pre- 
cisely localized on these emerging (quasiattractor) struc- 
tures, and ray-wave correspondence is fully established. 
We shall see below in Sec. IHII that it is the special struc- 
ture of the selection rule, namely the coupling of different 
maps, that induces the quasiattr acting features. 
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In this paper, we study the characteristics of quasiat- 
tractors emerging from simple coupled maps in order to 
gain a heuristic understanding of attracting phenomena. 
To this end, we introduce in Sec. [II] a toy model - a com- 
bination of standard and linear mappings - that clearly 
demonstrates the existence of the quasiattracting phe- 
nomenon and its origin. In Sec. IIIII we apply the insight 
gained to the optical system consisting of two dielectric 
disks, and explicitly show the mapping rules and their 
relation to the quasiattractor. We use the dependence 
of the selection rule on the refractive index of the disks 
in order to illustrate the parametric dependence of the 
emerging structures. Finally, we summarize our results 
in Sec. ED 



II. ORIGIN OF QUASIATTRACTING 
PHENOMENA 

A. Case study: Combined standard-linear mapping 

In order to understand the appearance of quasiattrac- 
tors, we introduce a map that combines the well-known 
standard and linear maps as an example of a piecewise 
smooth area-preserving map discussed in Refs. [7] as ori- 
gin of quasidissipative behaviour. The standard map 
is, physically speaking, a discrete-time analogue of the 
equation of the vertical pendulum and a prototype two- 
dimensional map commonly used in the study of various 
nonlinear phenomena of conservative systems P, LLOj . It 
is given by 






= In 



J n+ i (mod 27r), 
K s'm6 n (mod 27r), 



(1) 
for (0 n ,J n ) £ A. 



where 6 is the angle of rotation, / is the conjugate mo- 
mentum, if is a positive parameter which determines the 
dynamics of the map, and A is a subset of phase space 
introduced in connection with the selection rule. Equa- 
tion (pQ) is the usual standard map if we exclude the selec- 
tion rule, i.e., apply the map to all (0 ni I n ). For K = 0, 
the momentum is constant and the angle increases lin- 
early. As K increases, the phase space becomes increas- 
ingly chaotic as can be seen in the Poincare surface of 
sections (PSOSs) in Fig. [Q For if = 4.0, cf. Fig. [[] (a), 
one (split) island of stability remains, whereas all islands 
have disappeared for K = 9.0, see Fig. [1] (d). 

In addition to the standard map we introduce a lin- 
ear map with parameters a, j3 that has to be determined 
when (0 n ,I n ) e A, 



On+i = On + ol (mod 2tt) , 
In+i = I n + P (mod 2tt), 



(2) 



for (0 n ,7n) E A. 



B. Origin of the quasiattractor 

Each dynamics of standard maps and linear maps 
have been well understood for many years, but the dy- 
namics of the combined map is not the case. The re- 
gion 0.0 < n < 2.0 and 0.0 < I n < 2.5 is chosen 
for A, indicated by the red box in Fig. [I] (b). With 
(a, j3) = (2.5,0.0) the region A is mapped into the re- 
gion B represented by the dashed box in Fig. DJb) Ym ^ ne 
linear map (2). After some transient time, for K = 4.0 
the initial distribution uniformly prepared over the whole 
phase space converges into the island structure as shown 
in Fig. [2(b). This behavior can be easily understood as 
follows. The points starting from the regions outside the 
converging islands eventually reach the regions A due to 
ergodicity of chaotic dynamics. Some points of A are 
mapped into the islands according to the linear mapping 
of (2) since the region B contains the islands. Once the 
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FIG. 1: (color online) PSOSs of standard mapping for (a) 
K — 4.0 and (d) K — 9.0. Different color sets represent 
those from different initial points. The PSOS of the combined 
standard-linear map; (b, c) for K — 4.0 and (e, f) for K — 
9.0. The region A is denoted as the red solid box in (b). 
The region B, denoted as the dashed green box, is chosen 
as (b, e) (a,p) = (2.5,0.0) and (c) (a, f3) = (4.0,2.5). (f) 
Forward iterated sets starting from the region A using only 
the standard map. Red, green, and blue sets represent the 
first, second, and third iterated sets, respectively. 



points are put on the converging islands, they cannot es- 
cape from the islands since it forms an invariant set. If 
they are not put on the islands, they chaotically wander 
again until they are mapped into the islands. Therefore, 
the whole phase space converges to the islands as time 
goes on so that the islands form a quasiattractor [7[. 

The structure of the (quasi) attractor depends on where 
the region A and B are located. When B is chosen as the 
dashed box shown in Fig.l (c) with (a, /3) = (4.0,2.5), 
the quasiattractor consists of the chaotic components as 
well as the islands. The chaotic part resembles a usual 
strange attractor of dissipative systems [H[ . Its structure 
will be discussed in detail below. 

For K = 9.0, no stable islands exist in the standard 
map. After some transient time, a fractal quasiattractor 
emerges for (a, j3) = (2.5,0.0) as shown in Fig. [He). In 
order to understand the detailed structure of the fractal 
quasiattractor, let us consider the evolution of the region 
A only by the standard map (1). The first, second, and 
third iterations of the region A are denoted as red, green, 
and blue areas in Fig. [TJf ) , respectively, which shows a 
typical stretching and folding structure of chaotic dy- 
namics. Due to the linear map (2), the region A cannot 
follow such an evolution but is compulsorily mapped into 
B by (2). It means that the final attracting structure is 
likely to lack the areas generated by chaotic evolution of 
the region A. In fact, the empty area of Fig. [TJe) looks 
quite similar to the collection of the first, the second, and 
the third iterations of the region A shown in Fig. DJf). 
Note that this argument is only true for a short time 
evolution. However, the dissipative nature of our system 



guarantees that the short time dynamics is enough to un- 
derstand its main feature. This mechanism can also be 
applied to Fig. QJc) although the island part is dealt with 
separately. 

The fractal quasiattractors shown in Fig. 1(c) and (e) 
form an invariant set with a non-zero measure, which is 
called a fat fractal. As a parameter such as the location 
of the region B varies, the fractal quasiattr actor loses 
its stability and at the same time the island suddenly 
becomes an attr actor. This is a typical characteristic of 
the so-called crisis bifurcation [12, LL3j . It is known that 
forming a quasiattractor the phase space decreases on 
average linearly in time rather than exponentially [14j . 



III. QUASIATTRACTOR IN COUPLED 
DIELECTRIC CAVITIES 



The coupled dielectric cavities, schematically shown in 
Fig. 2, have a prominent feature [8|; when the ray es- 
capes from the cavity, it may either go to the infinity 
or collide with the other cavity so as to enter it. For 
the former case we force the ray to be reflected back 
into the original cavity rather than allowing the ray to 
escape from the system. For the latter case there ex- 
ists a finite probability that the ray can also be reflected 
from the surface of the other cavity. However, the ray is 
then forced to enter the cavity according to the rules set 
up in the RMDS. These forced operations play a role of 
the linear map (2) in the previous section. This shows 
how the RMDS works. Although this model is simple, 
it notably possesses the main feature of a quasiattrac- 
tor that we introduced in the previous section: The rays 
forcibly reflected back into the cavity, otherwise escap- 
ing from the cavity, indeed form quasiattr acting islands. 
Moreover they are shown to be directly associated with 
resonant modes obtained from wave calculation. 

In this section, we provide how the RMDS of two cou- 
pled dielectric disks is related to the coupled standard- 
linear map introduced in the previous section. We mostly 
exploit the results obtained from the scattering problem 
of two disk hard wall billiard outside the cavity [15] with 
Snell's law to incorporate the transmission across the cav- 
ity boundaries. 



A. System, selection rule, and mappings 

Two coupled disks schematically shown in Fig. 2 are 
described by geometric parameters, Rr, Rl, and D, 
which denote the radius of the right and the left disk, 
and the inter-disk distance, respectively, and by mate- 
rial parameters, namely the index of refraction of the 
right (tir) and the left (til) disks. For simplicity we set 
wr/ no = ^l/tio = n, with the index of refraction of 
the medium outside the cavities, no- The ray dynamics 
is traced by using Birkoff coordinate of (4>r(l), Or(l))i 
where (/)r(l) an d 6r(l) denote the azimuthal angle and 




FIG. 2: (color online) Schematic picture of the coupled two 
disks. The red arrows represent a typical trajectory of a ray 
transmitted from the left (L) to the right (R) disk. Note that 
in the left disk <j>l and Ol increases in the counterclockwise 
direction, while in the right the opposite is true. For example, 
in this figure, </>£,, <j>r, Or, and Or are positive and Ol and 0\ 
are negative. Rl and Rr is the radius of the left and the right 
disk, respectively, and D is the distance between two disks. 
til, riR and no are the indices of refraction of the left disk, 
the right disk and the outside, respectively. 



the angle of incidence in the right (left) disk, respectively. 
Notice that the sign of these angles are defined differently 
(See Fig. 2). 

The deterministic selection rule [8] is given as follows. 
For the sake of convenience, we assume the two disks are 
identical, i.e., R = Rr = Rl which allows us to find an 
analytic formulation of the problem. We will also con- 
sider what happens when two disks are not identical be- 
low. First of all, if the incident angle is larger than the 
critical angle C = arcsin(l/n), rays are totally reflected 
from the boundary and thus cannot escape from the disk. 
Transmission of rays from one disk to the other is possi- 
ble only when |0| < C . In addition, the escaped ray can 
arrive at the other disk only if \(j)\ < tt/2 is satisfied. As 
a matter of fact, the precise condition that the ray can 
be transmitted from one disk to the other is obtained by 
considering a simple geometry [15j and the Snell's law as 
follows 

arcsin(— 1/n) < < —9 m in for <p c < <p < tt/2, (3) 

-tax < < -6> min for - <j> c < (j)< </) c , 



< < arcsin(l/n) for - ir/2 < <j> < 



where 



arccos 



2R 



2R + DJ ' 

#min = arcsin (sin(6^ in )/n) , 
0max = arcsin (sin(<9^ ax )/n) , 



(4) 

(5) 
(6) 



with 



9 T 



9 T 

max 



/ Rsm</> \ . (R\ 

arcsin — ; ) — arcsin I — 1 , (7) 



arcsin 



I 

Rsim 



arcsin ( — ) , (8) 



where I = ^{2R + D) 2 + R 2 - 2R(2R + D) cos 0. The 
area of phase space satisfying the condition (3) is marked 




FIG. 3: (color online) The PSOSs of the two-disk map (TDM) 
for (a) n = 3.3, (b) n = 2.0, (c) n = 1.8, and (d) n = 1.5 
with D = 0.1 and Rr = Rl = 1.0. The right insets of (a)- 
(c) represent the corresponding enlarged PSOS of the part 
of the islands designated by the green boxes. The left insets 
represent the trajectories in real space corresponding to (a)- 
(c) the torus denoted as the thick red curve in the islands 
and (d) the chaotic quasiattractor. The outgoing (also called 
the region C) and the incoming area are represented as brown 
(dark gray) and yellow (light gray) area, respectively. They 
are symmetric over = and — due to the symmetric 
geometry shown in Fig. 2. 



by brown color in Fig. 3 and referred to as the outgoing 
area (also called the region C). The rays located in the 
outgoing area of one cavity are directly mapped into the 
corresponding area of phase space of the other cavity, 
which is called the incoming area marked by yellow color 
in Fig. 3. This constitutes the heart of the selection rule 



playing an equivalent role of the linear map (2) in that 
the outgoing and the incoming area directly correspond 
to the region A and B of the linear map (2), respectively. 
The incoming area should form mirror image of the out- 
going area due to the time-reversal symmetry and the 
reflection symmetry of the system. 

For completeness, we provide the full details of the 
map. Let us consider first the case that the rays stay 
within one of the disks. This is called the inside-disk 
dynamics, which is determined by the following simple 
mapping, 



9n+l = n . 



(9) 



When the ray enters the outgoing area of one disk by 
satisfying the condition (3), it is transmitted to reach the 
other disk. This transmission of the ray from one cavity 
to the other is described as the mapping 

</Wi = f{<i>n,0n) = T ((/> n ,arcsin(nsm(/9 n ))), (10) 

= arcsin(sm(T#(0 n ,arcsin(nsin(6VJ)))/n) , 
where T^ and Tq are given by |15j 



bn,0 n 



9n+l — 

= arcsm I sin 

Cl = Wn,C) 

'2R + D 



(11) 



sin (#n + <t>n) ) , 



= arcsm 



R 



R 



sin(C + 0n) - sin(C) , 



with 



\((j),0 T ) = (2R + D)cos(6 7 



Rcos0 T 



(12) 



[(2R + D) cos(<9 T 



Rcos6 T ) 2 - (2R + Df + 2R(2R + D) cos0. 
I 



The map of RMDS is thus summarized: if the condition 
(3) is satisfied Eq. ([TO]) is applied, otherwise Eq. (J9j). We 
call it two-disk map (TDM). 



B. Quasiattractor in coupled dielectric disks 

We obtain PSOS of the TDM after some transient 
time for various n and R with D = 0.1 using the uni- 
formly distributed initial points (</>o,#o) m the open re- 
gion, — arcsin(l/n) < < arcsin(l/n), of phase space. 
Figure [3] shows the PSOSs for several n with D = 0.1 
and Rr = Rl = 1.0, where the rays converge into the is- 



lands [Fig. 3(a)-(c)] or the fractal structure [Fig. 3(d)] so 
as to form the quasiattr actors. They look similar to those 
found in the standard-linear map of Sec. II. Figure 4(a) 
and (b) show the patterns of two resonant modes for 
n = 2.0 and n = 1.8, respectively, which are obtained nu- 
merically using the boundary element method [16|. The 
patterns are localized on the islands of Fig. 3(b) and 
(c), and closely resemble the ray trajectories presented 
in their insets. It clearly shows that the resonant modes 
are localized on the quasiattr actors. 

It is noted that the PSOS obtained here, Fig. 3, is 
slightly different from the usual PSOS in that the points 
are selected once the ray collides with the boundary ir- 




FIG. 4: Patterns of resonant modes obtained from wave cal- 
culation when (a) n — 2.0 and (b) n — 1.8, which are directly 
associated with the islands of Fig. 3(b) and (c), respectively. 
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FIG. 6: The PSOSs of the asymmetric TDM of Rr = 1.1R L 
obtained in the left disk, (a) n = 3.3, (b) n = 2.0, (c) n = 1.8, 
and (d) n = 1.5 with D = 0.1. 
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FIG. 5: (color online) The first (red), second (green), third 
(blue), and fourth (brown) iterated sets according to Eq. © 
starting from (a) the outgoing (black) and (b) the incoming 
area (black). 



respective of whether it is incident to or emerging from 
the boundary. Usually the point is chosen only when the 
ray is either incident or emerging. This explains why two 
islands are overlapped to form V shape in Fig. 3(a). We 
will discuss this point below. As n decreases the islands 
of the quasiattractors experience bifurcation; Firstly the 
island corresponding to the horizontal bouncing ball-like 
periodic orbit shown in the left insets of Fig. 3(a) and 
(b) becomes unstable so as to be split into the hexagonal 
shape periodic orbit shown in the left inset of Fig. 3(c). 
Secondly the hexagonal periodic orbit also loses its stabil- 
ity to form the fractal quasiattractor as show in Fig. 3(d). 
The islands in Fig. 3(a)-(c) consist of three parts, 
namely the parts in the outgoing (R ), the incoming (Ri), 
and the reflection region (R r ) in which all the ray is to- 
tally reflected in around \(j)\ = tt. The rays starting from 
R r on the island quasiattractor is mapped to those in R 01 
and consequently to Ri. This process is repeated, that 



is, R r — » R — » Ri — » R r — » R — » Ri. If time goes back- 
ward, outgoing and incoming rays are exchanged but the 
ordering rule is not affected. This means that the ray dy- 
namics on the islands is reversible in time, in agreement 
with the symmetric structures of the islands. However, 
in the case of a fractal-chaotic quasiattractor of Fig. [3] 
(d), the time reversal symmetry is broken since in gen- 
eral there can exist two possible origins for the ray inside 
the cavity; that is, one from the same cavity by reflec- 
tion from the boundary and the other coming from the 
other cavity by refraction. As shown in the standard- 
linear map of the previous section, here the fractal quasi- 
attractor also contains empty forbidden regions related 
to the short time successive iterations of outgoing and 
incoming regions. Figures [5] (a) and (b) show the first, 
second, third, and fourth forward iterated sets according 
to Eq. (J9]) from the outgoing and incoming areas, respec- 
tively. Iterated sets of the outgoing area in Fig. [5] (a) pro- 
duce the white forbidden region in Fig. [3] (d) whereas the 
iterations of the incoming area clearly leave their trace 
on the fractal quasiattractor. 

So far we have considered the identical disks, i.e. 
Rr = Rl, but our results can also be applied to the 
case of Rr ^ Rl- Figure [6] shows the PSOSs for various 
n with Rr = I.IRl, in which two disks are no longer 
identical. The PSOS of the left thus differs from that 
of the right. In Fig. 6 the PSOS is taken from the left. 
D = 0.1 is chosen as before. The PSOS is calculated by 
ray tracing method considering RMDS since the analytic 
expressions of TDM obtained above is only applicable for 
the identical disks. For n = 1.8 the quasiattractor of the 
hexagonal-shape periodic orbit is observed at Rr = Rl 
[Fig. 3(c)], but disappears at Rr = I.IRl and the frac- 
tal quasiattractor takes place [Fig. 6(c)]. If D varies, the 
topological structure of the PSOSs and the quasiattrac- 
tor are also changed (not shown here). The variations 
of the structures in the PSOS depending on system pa- 
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FIG. 7: The generalized PSOSs constructed by considering 
the (a, e) incident ray inside, (b, f) emerging ray inside, (c, g) 
incident ray outside, and (d, h) emerging ray outside cavities, 
which is schematically represented in each inset in the upper 
row. We choose n = 2.0 (the upper row) and n — 1.5 (the 
lower row) with D — 0.1 and Rr — Rl — 1-0. 



These generalized PSOSs are also useful to study ray- 
wave correspondence in coupled dielectric cavities since 
the corresponding quasi-eigenmode obtained from wave 
calculation can be represented by the generalized Husimi 
functions which also have four possible realizations in the 
exactly same manner [211 ]. 



IV. SUMMARY 

We have studied a quasiattracting phenomenon in cou- 
pled dielectric cavities. Its origin is in the ray-splitting 
dynamics modelled as the map based upon RMDS. The 
key feature of this map can be understood by considering 
a simple toy model, the standard-linear map. 
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The Birkoff coordinate of PSOS of dielectric cavities 
can be chosen as one of four possible components; either 
the incident or the emerging ray either inside or outside 
cavity. As far as a single cavity is concerned, all com- 
ponents are basically equal to each other since they are 
directly interconnected by the law of reflection and the 
Snell's law [191 l20j. However, in coupled dielectric cav- 
ities, every four components must be separately taken 
into account because of more complex ray dynamics that 
includes refraction and subsequent reentry into another 
component as well as reflection at the dielectric boundary 
back into the same component. 

We therefore introduce a generalization of the PSOS 
that comprises incident and emerging rays inside and 
outside cavities. Figure [71 (a)-(d) show the generalized 
PSOSs of our ray model for n = 2.0, D = 0.1, and 
Rr = Rl = 1.0. As mentioned above, the PSOS of 
Fig. 3 was constructed such that the points are selected 
once the ray collides with the boundary. It means that 
Fig. 3(b) is the combination of Fig. 7(a) and (b), which 
is the PSOS chosen from the incident and emerging ray 
inside the cavity, respectively. Considering the nature of 
the coupled two disks, whose ray is either reflected from 
or transmitted through the boundary, the appropriate 
PSOS should contain both cases. This is the reason why 
the PSOSs was made in this way. The case of chaotic 
PSOSs with fractal quasiattracting structures is shown 
for n = 1.5, D = 0.1, and Rr = R L = 1-0 in Fig. (e)- 
(h). Again, Fig. 3(d) is obtained from merging Fig. 7(e) 
and (f). 



Appendix: Monodromy matrix of hexagonal-shaped 
periodic orbit 

The hexagonal-shaped periodic orbit (HSPO) is de- 
scribed by 12 successive mappings consisting of 6 trans- 
lations (4 inside the cavity and 2 outside), 4 refractions 
across the boundary and 2 reflections at <p = tt. The mon- 
odromy matrix of the HSPO can then be constructed as 
multiplication of the monodromy matrices of each map- 
ping. Note that the monodromy matrix of 2 reflections is 
an identity so that it can be ignored. For R = Rr = Rl, 
the incident angle of the hexagonal periodic orbit is given 
by 6 = arccos(n/2) when n < 2.0. There is no hexagonal- 
shaped periodic orbit if n > 2.0. The monodromy matri- 
ces of the translation inside and outside the cavity are 



Mi 



2R 
cos 
1 



and 



M = 



, Rcos 20 
I 2 cos 20 



.Rcos 20 



(13) 



(14) 



respectively, I = D + 2R(1 — cos 2$). 

The monodromy matrices of the refraction from inside 
to outside and the opposite can be obtained from Snell's 
law. 



Mi 



1 
n 



(15) 



and 



Mb 2 



1 
- 1 



(16) 



respectively. 

Finally, the monodromy matrix of the hexagonal- 
shaped periodic orbit is 

M = MiMb^oMb^iMiMb^qMb^i. (17) 



The stability of the hexagonal-shaped periodic orbit de- 
pends on the value of Tr M. For R = 1.0 and D = 0.1, 
in the case of Tr M > 2 when n < 1.688, the periodic 
orbit is linearly unstable but in the case of Tr M < 2 
when n > 1.688, the periodic orbit is linearly stable. 
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